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ABSTRACT 

The X-ray luminosity and temperature of clusters and groups of galaxies do not scale in a self-similar 
manner. This has often been interpreted as a sign that the intracluster medium has been substantially 
heated by non-gravitational sources. In this paper, we propose a simple model which, instead, uses the 
properties of galaxy formation to explain the observations. Drawing on available observations, we show 
that there is evidence that the efficiency of galaxy formation was higher in groups than in clusters. If 
confirmed, this would deplete the low-entropy gas in groups, increase their central entropy and decrease 
their X-ray luminosity. A simple, empirical, hydrostatic model appears to match both the luminosity- 
temperature relation of clusters and properties of their internal structure as well. 

Subject headings: cosmology: theory, intergalactic medium 



1. INTRODUCTION 

Clusters and groups of galaxies are composed of galax- 
ies, hot X-ray emitting gas, and a gravitationally dominant 
dark halo. Although this basic picture is well-understood, 
there remain a number of puzzles that prevent clusters 
from being fully used as standard candles. For example, 
simple scaling relations (and detailed simulations) predict 
that the bolometric X-ray luminosity should scale with the 
temperature of the gas according to Lx ~ T 2 , while obser- 
vations indicate Lx ~ T 3 (Kaiser 1991; Evrard & Henry 
1991; Navarro, Frenk & White 1995; Bryan & Norman 
1998). Early on it was suggested that groups might have 
a lower baryon fraction because of more efficient star for- 
mation (e.g. David & Blumcnthal 1992; Thomas & Couch- 
man 1992; Bower 1997), but later thinking has focussed 
on the idea of additional heating of the gas, often assumed 
to be from supernovae or active galactic nuclei (e.g. Kaiser 
1991; Bower et al. 2000). This non-gravitational heating 
would decrease the central density and, because the X-ray 
emissivity is proportional to the density squared, reduce 
the luminosity. Because of the lower pressures in smaller 
clusters and groups, this would preferentially affect them, 
steepening the Lx — T relation (Cavaliere, Menci & Tozzi 
1997). This viewpoint was strengthened by the discovery 
of an apparent entropy "floor" in the centers of groups and 
clusters (Ponman, Cannon & Navarro 1999). This is con- 
sistent with the idea of a source of heat which raises the 
entropy of the gas to a fixed, minimum level; smaller clus- 
ters have lower entropies and so are more affected than 
large clusters. This model has been developed in some 
detail in a number of papers (Balogh et al. 1999; Loewen- 
stein 2000; Wu, Fabian & Nulsen 1999; Cavaliere, Menci 
& Tozzi 1999; Valageas & Silk 1999) and appears to be 
capable of naturally explaining the observations. 

The amount of heating required is substantial. Although 
estimates vary, it seems likely that about 1 keV per parti- 
cle is needed (Lloyd-Davies et al. 2000), an amount which 
may be challenging to explain from supernova heating 
alone (Kravtsov & Yepes 2000). Another difficulty is that 
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observations of the Lya forest indicate a much lower tem- 
perature for the majority of the intergalactic medium at 
z ~ 2 - 3 (Bryan & Machacek 2000; Schaye et al. 1999), 
a condition which may extend to even lower redshift (Ri- 
cotti, Gnedin & Shull 2000). Although hardly conclusive, 
these concerns may be pointing toward another explana- 
tion for the observations. 

This paper argues that cooling and the resultant galaxy 
formation are sufficient, by themselves, to explain all of 
these observations and that substantial heating is not re- 
quired. We draw mostly on two simple ideas: (1) small 
clusters and groups have converted more of their baryonic 
gas into galaxies than have large clusters and (2) the gas 
which goes to form the galaxies is preferentially lower in 
entropy, thus raising the mean entropy of the gas which 
remains. This means that not only do small clusters have 
a smaller gas fraction {f gas ), the gas which is there has a 
higher entropy — and lower density — than it would in 
simple self-similar scaling models. Because of the density- 
squared nature of the X-ray emission, this substantially 
diminishes the luminosity of groups and small clusters, re- 
sulting in a steeper Lx — T relation, as observed. The 
effective entropy increase is most noticeable in the cen- 
ter of the cluster, which is just where the entropy floor 
is observed. In what follows, we explore a simple model 
to investigate if this hypothesis can match the observa- 
tions quantitatively. We will also show that there is some 
empirical support for the first assumption. 

2. THE MODEL 

The model described here is built on the assumption 
that galaxy formation is not uniformly efficient in all envi- 
ronments. Since theoretical arguments can and have been 
made both ways, we attempt to address this point with 
observations drawn from the literature. We searched the 
literature and found three studies that computed stellar, 
gas and total masses within the same radius (many more 
computed a subset of these three quantities but we only 
used those that computed all of the quantities to insure 
self-consistency). Mulchaey et al. (1996) used their own 



ROSAT and optical observations combined with other re- 
sults from the literature to compile a list of 16 groups 
with masses computed out to Rx, the maximum radius 
at which X-ray emission could be observed. Hwang et 
al. (1999) used ASCA observations of five intermediate 
mass systems and also computed masses out to Rx- Fi- 
nally, Cirimele et al. (1997) studied 12 Abell clusters with 
ROSAT, and tabulate masses computed out to 1.5 Mpc 
(which is close to Rx for their clusters). All studies used 
similar, although not identical, stellar mass to light ratios 
(corrected for morphological variations) and all used the 
same cosmological parameters. 

In Figure 1, we plot the relative stellar and gas masses 
from these three studies. This shows that the hot gas com- 
ponent dominates over galaxies in the most massive clus- 
ters of galaxies. For smaller systems, the scatter increases 
significantly; however, there is a trend towards an increas- 
ing stellar contribution and decreasing gas contribution 
for lower mass clusters and groups. To make this clearer, 
we divide the sample into four equal-sized groups, ordered 
by temperature, and plot the median for each group. To 
check the statistical strength of the trend, we divided the 
sample into two groups (those with temperatures below 
and above 2 keV), and separately computed the median 
galaxy mass ratios. Then, using median statistics (e.g. 
Gott et al. 2000), we find the probability that the median 
of the high temperature group is smaller than that of the 
low temperature group to be 0.986. The gas mass ratio 
trend was even more significant. 
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Fig. 1. — The ratio of mass in hot gas (solid symbols) and in 
stellar systems (open symbols) to total gravitating mass as a func- 
tion of temperature. The circles are from Cirimele, Nesci & Trevese 
(1997), the squares from Hwang et al. (1999) and the diamonds from 
Mulchacy et al. (1996). Solid and dashed lines arc median values 
of the binned distribution for the gas ratio and stellar ratio, respec- 
tively. The dot-dashed line is the model discussed in the text. All 
observations have been adjusted to h = 0.65. 

The most straightforward explanation of the trend in 
Figure 1 is that the efficiency of galaxy formation varies 
from groups to clusters. This is the hypothesis that we will 
examine in this paper, but there are certainly other expla- 
nations. For example, it is possible that the relative mix 



of gas and stars changes outside of the measured region 
(i.e. R > Rx), which is generally a smaller fraction of the 
virial radius for groups than for clusters. Moreover, it is 
difficult to be conclusive for a heterogeneous sample of this 
sort, since the groups and clusters were examined by dif- 
ferent authors using slightly different methods. However, 
all of the studies did use the same basic methodology, and 
adopted similar parameters. Also, the trend itself does not 
depend on correctly determining the total mass, since the 
ratio M gas /M op t — which is independent of total mass — 
also decreases with temperature. 

There are a number of other pieces of evidence which 
support this basic conclusion. For example, Arnaud & 
Evrard (1999) and Mohr, Matheison & Evrard (1999) find 
a trend of decreasing hot gas fraction with decreasing tem- 
perature for their samples. This pattern has sometimes 
been taken to imply that gas has been ejected from smaller 
clusters and groups, despite the large amount of energy 
required to do this. However, it seems equally possible 
that this gas is in the form of stars. This would also 
be consistent with a higher mass-to-light ratio for large 
clusters than for groups (Girardi et al. 2000; Adami et 
al. 1998; Hradecky et al. 2000; Ramella, Pisani & Geller 
1997), although see David, Jones & Forman (1995). Weak 
lensing of groups should provide useful constraints; pre- 
liminary results indicate the the mass-to-light ratios of 
groups are somewhat lower than clusters (Hoekstra et al. 
2000). Finally, constraints from galaxy clustering indicates 
the number of galaxies in a halo must grow more slowly 
then the mass of the halo (Scoccimarro et al. 2000; Scljak 
2000), consistent with the trend presented here. Despite 
this circumstantial evidence, we cannot prove that the effi- 
ciency of galaxy formation depends on environment; all we 
can do is show that the available data are consistent with 
the trend shown in Figure 1. In the rest of the paper, we 
will assume this to be true and examine the consequences 
that follow. 

In Figure 1, we show the relation, 
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which we will take to be the stellar mass fraction in this 
paper, and is the result of minimizing the absolute de- 
viation of the mass ratio. The gas fraction is simply 

fgas = fbaryon ~ fstar, where fbaryon = 0.16 is Compatible 

with the cosmological model we have chosen. Specifically, 
this is a flat model with (D, ,Cl b ,h) = (0.35,0.056,0.65), 
where f^o is the ratio of the mass density to the criti- 
cal density and h is the Hubble constant in units of 100 
km/s/Mpc. The results are most sensitive to the value of 
the Hubble constant since the ratio of stellar to gas mass 
varies as /i -3 / 2 ; the other parameters play almost no role. 
In order to build a concrete model for the structure of 
a group or cluster of galaxies we assume that: (1) the 
clusters are spherically symmetric and in hydrostatic equi- 
librium; (2) the hot gas and stellar fractions are as given 
above; (3) the gas which is converted into stars comes from 
the lowest entropy gas in the cluster and all other fluid el- 
ements lie on the same adiabat they would have without 
cooling or star formation. From experience gained with 
numerical simulations, we know that while clusters are not 
in exact hydrostatic equilibrium this assumption is a rea- 
sonable approximation. The second assumption has some 
empirical basis, as previously discussed. We will return to 



a discussion of the last assumption. 

To create a cluster of a given mass M, we assume the 
dark matter density is described by (Navarro, Frcnk & 
White 1996): 

(2) 



p(x) 200 



po 3 ln(l + c) -c/(l + c) .t(1 + cx) 2 

where po = 3H 2 /8irG is the critical density, a; = 
r/r2oo and c is a concentration parameter that depends 
weakly on mass in the range of interest. We take c = 
8.5(Mft/lO 15 M )- 0086 . The radius r 200 is defined by 
M = 800p 7rr| 00 /3. 




Fig. 2. — The gas profiles of clusters after galaxy formation for 
temperature ranging from 10 kcV (dotted) to 1 keV (dot-dashed). 
The top (solid) curve shows the initial profile. Note that the decrease 
with cluster temperature is entirely due to the increased efficiency 
of galaxy formation — without this all curves would look like the 
initial profile. 

The gas distribution without galaxy formation is as- 
sumed to have the same distribution as eq.(2). The tem- 
perature profile is determined by solving the equation of 
hydrostatic equilibrium for a spherical profile: dP/dr = 
—pbGM(r)/r 2 . We assume that the gas does not con- 
tribute to the gravitational potential and adopt a pressure- 
free external boundary condition (i.e. P = at s = 1). 
The result matches the density and temperature distri- 
bution in numerical simulations over the vast majority of 
the cluster (e.g. Frenk et al. 2000). It fails in the very 
center (where the numerical models are the least certain) 
but this represents a small fraction of the mass (< 1%) 
which will end up being converted to galaxies anyway. An 
isothermal temperature profile produces results which are 
broadly similar; however, this density distribution does not 
well describe the simulation results, particularly at large 
radius (x > 0.3) where much of the mass resides. 

Once the no-cooling cluster has been constructed, we 
can then compute the structure of the cluster including 
galaxy formation. The equation of hydrostatic equilibrium 
sets the pressure distribution, but we need one more con- 
straint to uniquely fix the density and temperature pro- 
files. This comes from the entropy (S = \nT/p b ) dis- 
tribution of the gas, which is a monotonically increasing 



function of radius. Since by our earlier assumption, the 
gas which cools into galaxies comes from the lowest en- 
tropy part of the distribution, galaxy formation is equiva- 
lent to removing from the center an amount of gas equal 
to f s tarM. The remaining gas is then distributed over the 
whole cluster, under the assumption that it docs not cool 
at all. The known entropy distribution of this gas com- 
bined with the equation of hydrostatic equilibrium is suffi- 
cient to specify the gas and temperature profiles uniquely. 
More precisely, we guess a central pressure and then inte- 
grate the equation of hydrostatic equilibrium outwards in 
radius, or cquivalcntly, enclosed baryonic mass Mb(r). At 
each point, the density and temperature is computed from 
the pressure and entropy S(Mb), with the entropy coming 
from the no-cooling case. We iterate this procedure with a 
new central pressure each time until a profile is produced 
which conserves mass (i.e. M(,(r2oo) = fbaryonM). Al- 
though this may superficially appear to be a cooling-flow 
model, it is not. The gas that cools into galaxies does so at 
high redshift, before the cluster even forms. The gas which 
goes into the galaxies is at the center of the smaller halos 
which merge to form a cluster. If that gas hadn't cooled 
to form galaxies, it would have ended up in the center of 
the cluster, since low entropy gas will sink to the center. 

Figure 2 shows the resulting density profiles for a range 
of cluster temperatures. The temperature profiles remain 
mostly flat, although they increase slowly towards the cen- 
ter and have a somewhat larger mean compared to their 
non-cooling equivalents. 
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Fig. 3. — The data points show the observed bolomctric X-ray 
luminosity temperature relation, with open points from Arnaud & 
Evrard (1998) and filled points from Hclsdon & Ponman (2000). 
The solid line is the Lx — T relation for the cluster model discussed 
in the text, while the dashed line is for a self-similar model. The dot- 
dashed curve shows the result of adopting the variable f ga s fraction, 
but using a self-similar profile for the gas. 



3. COMPARISON TO OBSERVATIONS 

The first test is to see if the resulting clusters can match 
the luminosity-temperature relation. For each cluster the 
X-ray luminosity and luminosity-weighted temperature is 



computed with a Raymond-Smith code (1977; 1992 ver- 
sion) assuming a constant metallicity of one-third solar. 
The result is shown in Figure 3 along with the result for 
a self-similar model (this is actually flatter than Lx ~ T 2 
due to the increased importance of metal lines at low tem- 
peratures). In the same figure, we also show the Lx —T re- 
lation which would result from keeping a self-similar shape 
for the gas profile, but reducing the gas fraction as speci- 
fied by eq. (1). 

The agreement is good, except at very low temperatures, 
where the observations fall below the theoretical curve. 
The surface brightness of these poor groups is very low, 
and so a significant fraction of their luminosity could be 
lost. Helsdon & Ponman (2000) estimate that for their 
lowest temperature cluster group (below 1 keV), the total 
flux is underestimated by about 40% relative to the higher 
temperature groups. Also, there is some evidence that 
the metallicity of small groups may be significantly less 
than that of larger groups (Davis, Mulchaey & Mushotzky 
1999). These two effects would reduce the predicted lumi- 
nosity by a factor of 2-3. 

There is another line of evidence which has been taken 
as strong evidence of preheating: Ponman, Cannon & 
Navarro (1999) show that the central entropy (at r = 
0-lr2oo) in clusters and groups does not scale in a self- 
similar fashion. Their data are reproduced in Figure 4, 
along with the curve predicted from the model described 
in this paper. The model matches the observed data. Also 
shown is the self-similar relation (with constant f ga s)- 
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Fig. 4. — The entropy of the gas at r = 0.1r200 for a collection 
of clusters and groups (Ponman, Cannon & Navarro (1999)). The 
solid line shows the prediction for the model described in this paper, 
while the dashed line is the self-similar relation. 

There is a range of other observations against which 
this model could be tested; we restrict ourselves to two 
others. There is some evidence that the profiles of groups 
is flatter than that of clusters. For example, if the X- 
ray surface brightness profile is fit with a beta model: 
S(R) = 5 (1 + (R/R c ) 2 )- 3/3+1 / 2 , then the outer slope 
/3 is ~ 0.7 for large clusters, ranging down to f3 ~ 0.4 
for groups (Mohr, Matheisen & Evrard 1999; Helsdon & 



Ponman 2000). In this expression, R c is the projected 
cluster core radius. As might be imagined from Figure 2, 
our model also shows this trend. After integrating along 
lines of sight, the resulting surface brightness profile is well 
fit by (3 = 0.75 for a 10.2 keV cluster and (3 = 0.5 for a 1.2 
keV group. 

The last check we make is to examine the evolution 
of the Lx — T relation with redshift. Observations 
(Mushotzky & Scharf 1997; Schindler 1999; Fairley et al. 
2000) show that there is very little change in this relation 
toz~ 0.5, although the amount of data are still limited. 
Unfortunately, no high redshift equivalent of Figure 1 ex- 
ists; however, if we assume that the ratios do not change 
appreciably then this models also predicts little evolution 
to modest redshifts. Indeed, most models which correctly 
predict Lx ~ T 3 will reproduce this lack of evolution. The 
reason is simple: for a fixed mass, the luminosity scales 
roughly as (1 + z) 3 as long as the profile doesn't change 
very much when expressed as a function of r/r2oo ■ Also for 
a fixed mass, the virial temperature scales as (1 + z) and 
so modifying z moves a cluster parallel to the Lx ~ T 3 
relation. 

4. DISCUSSION 

In this letter, we have described a simple model of clus- 
ter formation which reproduces the self-similar breaking 
observations without recourse to non-gravitational heat- 
ing. There are two key assumptions in this model. The 
first is that galaxy formation was more efficient in groups 
than in clusters; as discussed in section 2 there is some em- 
pirical evidence for this. Certainly the morphology-density 
relation shows that galaxies are sensitive to their environ- 
ment. From a theoretical standpoint, this could arising 
from biasing (David & Blumenthal 1992) or from cooling 
and shocking of gas. 

The second important assumption is that the lowest 
entropy gas is converted into galaxies, while the high- 
entropy gas retains the same entropy it would have had 
without galaxy formation. Clearly this is an idealization: 
in practice the rest of the gas will suffer some radiative 
losses (and if it cools substantially, this will invalidate the 
model assumed here). However, the approximation is self- 
consistent in that the remaining gas has a cooling time 
comparable to or longer than the Hubble time. It is also 
true that cooling of hot gas tends to occur catastrophi- 
cally (e.g. Thoul & Weinberg 1995). That is, it remains 
hot with little cooling until it passes through a cooling 
front, where the density suddenly increases by orders of 
magnitude. The gas at large radius moves towards the 
center without changing it's entropy. It seems clear that 
numerical simulations will be required to test these argu- 
ments, although it will be computationally challenging to 
do so. It is possible that some of the effects described 
in this paper may have already been seen in simulations 
(Pearce et al. 2000). 

We acknowledge useful discussions with Mark Voit, Joe 
Mohr and the referee, Richard Bower. Support for this 
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grant HF-01104.01-98A from the Space Telescope Science 
Institute, which is operated under NASA contract NAS6- 
26555. 
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